seis=seis*diag(scales);
